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Very long transients are found in the partially ordered phase of type II of globally coupled logistic maps. The transients always 
lead the system in this phase to a state with a few synchronous clusters. This transient behaviour is not significantly influenced 
by the introduction of weak noises. However, such noises generally favor cluster partitions with more stable periodic dynamics. 
PACS number(s): 05.45.-a, 05.45.Xt, 05.40. Ca 



Globally coupled maps (GCM) formed by ensembles of 
logistic maps have been used as a paradigm of complex 
collective dynamic behaviour for a decade Origi- 
nally, GCM were introduced as a mean-field approach to 
coupled map lattices || , but later the nontrivial dynam- 
ics and the rich collective phenomenology displayed by 
that system made it a subject worth of study in itself. 
One of the main properties of globally coupled logistic 
maps is the presence of different phases characterized by 
turbulent (non-synchronized) behaviour, clustering, and 
global synchronization . The formation of a number of 
subgroups of synchronized elements out of a symmetri- 
cal ensemble has a high relevance for many applications, 
such as the organization of the immune or the neural sys- 
tem, ecological networks, cell differentiation, and struc- 
turing of social hierarchies. Therefore, the GCM phases 
in which the system displays clustering have been inten- 
sively studied 

The former studies about the complex collective behav- 
ior displayed by GCM are usually based on the classical 
classification of the phase-space introduced by K. Kaneko 
As we report here, the partially ordered phase of type 
II is equivalent to its neighboring ordered phase, with 
the only difference that very long transients preceed the 
achievement of the final attractor. This implies a revi- 
sion of the phase space of GCM, and shows that there 
is a strong non-monotonous dependence of the transient 
length with the system parameters, which has to be taken 
into account in any numerical study. 

The simplest globally coupled discrete-time system is 
given by 

N 

Xi (t + 1) = (1 - e)f( Xi (t)) + ± £ /foW) • 

where the individual element evolves according to the lo- 
gistic map f(x) = 1 — ax 2 , N is the total number of maps 
and e specifies the coupling strength. In general, an at- 
tractor of this dynamical system is formed by a number 
K. of synchronous clusters each containing Nh elements, 
k = 1, . . .K, and can be characterized by means of the 
partition (/C; N% > N2 > . . . > N/c). For convenience, 
this classification also includes one-element "clusters" 
(iVfe = 1) that actually correspond to individual non- 
entrained elements. Thus, the partition (N; 1, 1, 1) 
corresponds to the asynchronous state of the entire en- 
semble, while the partition (1; N) represents its fully syn- 
chronous state, where all N elements belong to a single 



cluster. In addition to these two states, the system would 
generally also have other partitions where a certain num- 
ber of synchronous groups of elements with TV > K, > 1 
are present. The choise of an attractor with a particu- 
lar partition is determined by the initial conditions. The 
attractor corresponding to each initial condition is char- 
acterized by a certain number /C m of clusters after a tran- 
sient has elapsed. 

The phase space of the GCM (Q) has been described 
by K. Kaneko using the average cluster number K, = 

. . , M enumer- 
Four different 



Em=i^m/^i where the index m = 1, . 
ates the set of employed initial conditions, 
phases have been identified 0]: 



1. Coherent phase 

jectory (x^t) = Xjfr^Vi, j 
synchronous cluster (/C = 1). 



The elements follow the same tra- 
t), forming a single 



2. Ordered phase. Almost all basin volume is occupied 
by a few-cluster attractor (/C is small and does not 
grow with N). 

3. Partially ordered phase. Coexistence of many- 
cluster and few-cluster attractors {JC is large and 
grows with N). 



Turbulent phase. No 
elements (/C = N). 



synchronization among the 



The coexistence of many-cluster and few-cluster at- 
tractors has been observed by K. Kaneko in two different 
parameter intervals. One of them separated the ordered 
and the turbulent phases. Here the system is in the par- 
tially ordered phase of type I, also called the intermittent 
phase. The other interval lies between the regions occu- 
pied by the ordered and coherent phases. In this inter- 
val the system is in the partially ordered phase of type 
II (called the "glassy" phase in the initial publication 
[||). The typical parameter intervals are 1.56 < a < 1.80 
for e = 0.3 (partially ordered phase of type II) P,^3|, 
and 1.58 < a < 1.69 for e = 0.1 (intermittent phase) 

ma. 

To compute the asymptotic properties of a dynami- 
cal system, one has to ensure that the system has had 
enough time to approach its final state, i.e. that the 
dynamical attractor for the given initial conditions has 
been reached. Slow relaxation is indeed known for some 
dynamical systems (see, e.g., pT|). The properties of the 
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transients of GCM have not yet been sufficiently inves- 
tigated. The aim of the present Letter is to systemat- 
ically study the transient behaviour of GCM, described 
by equation (Q). Our principal result is that inside the 
whole parameter region, corresponding to partially or- 
dered ("glassy") phases of type II, only few-cluster at- 
tractors are observed after very long transients. This 
result holds also when weak noises are added. It contra- 
dicts to what has previously been reported by K. Kaneko 
& 

We performed long runs of up to T = 10 iterations 
and recorded the time at which the final minimum value 
of /C was reached in the parameter region corresponding 
to the partially ordered phase of type II. To this end, 
the partition (/C; Ni > JV-j > • • • > N/c) was determined 
every At = 25 — 50 time steps. Double precision real 
numbers were used in these computations, ensuring the 
absolute precision of 10 -16 . Two elements were taken to 
belong to the same cluster only if they had exactly the 
same state within the double computer precision, i.e. if 
\ Xi (t) - Xj (t) | < 10- 16 . 
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FIG. 1. Average time required to approach an attractor as 
function of the control parameter a for the coupling strength 
e = 0.3 in a system of size N = 100. Averages over 100 dif- 
ferent initial conditions are performed. Vertical bars indicate 
minimum and maximum values of the transient time. 

The computed average transient length for e = 0.3 as 
a function of a is shown in Figure 1. We see that the 
transients may extend up to tens of thousands and even 
millions of time steps. They become especially long near 
a = 1.6. Previous numerical studies (l],g,[l3| were limited 
to much shorter evolution times (up to 10 4 iterations) 
and therefore some of the behaviour observed in these 
studies essentially corresponded to transients. This be- 
comes clear if we compare our Fig. 1 with Fig. 9 in Ref. 

: A strong increase in the mean number of clusters K, 
was reported exactly where the transient length greatly 
increases (exceeding 3000 time steps). Our investigation 
reveals that, after long transients, only attractors with 



K < 2 are typically found for a < 1.65, and only at- 
tractors with JC < 6 are observed for 1.65 < a < 2. 
Similar results are also obtained in our calculations for 
e = 0.25, 0.35, and 0.4. 
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FIG. 2. Time dependence of the number of clusters in a 
system of size N = 100 for the coupling intensity e = 0.3 
at a — 1.6 (main plot) and a = 1.8 (inset), and 30 different 
initial conditions in each case. 

The time dependence of the number of clusters IC(t) 
at e = 0.3 for 30 different initial conditions is shown in 
Figure 2 for a — 1.6 (main plot) and a = 1.8 (inset). For 
a = 1.6, the number of clusters is indeed large during the 
initial evolution, and comparable with the total size of 
the system (N — 100). Later on the number of clusters 
is slowly decreasing and eventually only attractors with 
K, = 2 (but different partitions N\, N2) are found at 
this value of the parameter a. The system evolution at 
a = 1.8 is essentially similar, though it is characterized 
by a much faster convergence to the final states (note 
the difference by almost two orders of magnitude in the 
time scales in these two plots). Another difference is that 
for a = 1.8 the final states with various cluster numbers 
K, = 2, 3 and 4 are observed. By averaging over a large 
number of initial conditions, the time dependence for the 
relaxation of the mean cluster number IC(t) to its asymp- 
totic value /C has been obtained. Figure 3 shows in the 
logarithmic scale the time dependence of the quantity 
5fC = K,(t) — K, for e = 0.3 and a = 1.6. We clearly see 
that the relaxation is exponential, 8K, oc exp{— (3i). 

We have further analysed how the mean transient time 
t — P^ 1 depended on the system size N. The explored 
interval of system sizes was 2 4 < N < 2 12 ; we have used 
several values of a and fixed the coupling strength at 
e = 0.3. We did not find any strong variation of r with 
N , i.e. the order of magnitude of r did not depend on 
the system size. The transient length depicted in Fig. 1 
is characteristic for almost three decades of variation in 
the system size N. 
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FIG. 3. Exponential relaxation to the asymptotic state for 
the coupling intensity e — 0.3 and three different values of 
the control parameter a. The slope of the curves (from top 
to bottom) is (3 = 1.1 x 1(T 4 , 2.4 x 1(T 4 , and 5.4 x 1(T 4 , 
and errors are of order 10 -7 . The system size is N = 100. 
Averaging over M = 10 3 initial conditions is performed. 

The presence of very long transients indicates that 
the system may be sensitive to the application of noise. 
Indeed, for the so-called Milnor attractors even a tiny 
perturbation would suffice to destabilize the asymptotic 
state [Q. The existence of Milnor attractors has been 
discussed both for the partially ordered phase of type II 
fl3f and for the intermittent phase fl(i|j . To analyze the 
effect of weak random perturbations, we have modified 
equations d) by adding a noise term rj ri(t). We have 
chosen a small noise intensity r\ = 10~ 10 ; independent 
random numbers Ti(t) £ (—1, 1) are drawn anew from a 
uniform distribution for each element and at each time 
step. Noise prevents the spurious synchronization of ele- 
ments in the system: If the states of two maps i and j are 
equal (with computer precision) at time t, they will follow 
identical trajectories for all t' > t in a pure deterministic 
system. When noise is added, spurious attractors are not 
attained and only robust attractors should be detected. 

In the presence of noise the states of elements in a 
cluster cannot be identical. To define a cluster, we have 
to choose a certain finite precision 5 and say that ele- 
ments i and j belong to the same cluster at time t if 
\xi(t) — xj(t)\ < S (cf. the respective definition for the 
case of randomly coupled maps We have found that 
the application of weak noise does not qualitatively in- 
fluence the above-described evolution. Figure 4 shows in 
logarithmic scale the mean number IC(t) of clusters as 
function of time when weak noise is present (all other 
parameters are the same as in Fig. 3). The system 
still evolves towards final distributions characterized by 
a few large clusters. Typically, a slower convergence to 
the limit value K,~ 2 was observed when the noise was 
acting. Only in a very narrow domain 1.60 < a < 1.62 
did noise seem to prevent convergence to a few-cluster at- 



tractor. Note that this area coincides with the maximum 
transient length in the deterministic case, and is near the 
boundary where the single synchronous cluster becomes 
unstable (it is known that the unstabilization of the co- 
herent phase proceeds through a power-law divergence of 
the transient lenght §). 
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FIG. 4. Exponential relaxation in GCM added with noise 
of amplitude r\ = 10~ 10 (other parameters as in Fig. 3). The 
precision is S = rj 1,/2 . The slope of the curves (from top to 
bottom) is (3 = 8 x 10" 5 , 2 x 10" 4 , and 5 x 10" 4 . 

The dynamics corresponding to a particular cluster 
partition in our simulations was either periodic or (intrin- 
sically) chaotic. To detect intrinsically chaotic dynamics, 
local Lyapunov exponents were examined. After a fixed 
transient of length T = 10 7 we numerically calculated 
the local Lyapunov exponent 

T+t N 

X m (e,a) = - ^5>g|/'(^(t))| (2) 

corresponding to the trajectories Xj(t) of elements j for 
the given initial condition m and parameters e and a. The 
averaging time was always t = 10 4 . Positive exponents 
correspond to chaotic dynamics. The same procedure 
was used both in the presence and in absence of the noise. 

When noise is acting, it may, in principle, induce tran- 
sitions from one cluster partition to another. Our nu- 
merical investigations show, however, that such transi- 
tions actually take place in the presence of very weak 
noises only if the dynamics corresponding to a particular 
cluster partition is intrinsically chaotic. This observation 
leads us to a conjecture that Milnor attractors in GCM 
are, perhaps, only generated by cluster partitions with 
intrinsically chaotic dynamics. Cluster partitions with 
periodic attractors are stable against a finite amount of 
perturbation, while the system leaves with certainty a 
partition with a chaotic attractor in a finite time when 
noise is present (see also p~0|] ) . 
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of few- and many-cluster attractors is indeed observed. 
These results are important for the general classification 
of dynamical behaviour in GCM. 

The authors acknowledge the financial support from 
the Alexander- von-Humboldt Foundation (Germany). 
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FIG. 5. Fraction / of initial conditions leading to a chaotic 
attractor for GCM with and without noise. The parameters 
are e = 0.3 and N = 100 and we have averaged over M — 10 3 
initial conditions. 

Interestingly, the addition of noise favors the attaine- 
ment of periodic attractors. Figure 5 shows the fraction 
/ of initial conditions leading to a cluster partition with 
(intrinsically) chaotic dynamics with and without noise 
for the coupling intensity e = 0.3, as identified by means 
of (||) . We see that this fraction is strongly reduced in the 
presence of noise around the parameter value a = 1.61. 
To explain this, suppose that the system has approached 
a cluster partition with chaotic dynamics. Given that 
any such partition is destabilized even by weak noise, we 
expect that elements would spend only some time near 
this attractor, but then one of them would change its 
cluster affiliation and a new cluster partition would thus 
be produced. As long as this new partition is also chaotic, 
the system again easily escapes and the same procedure 
repeats until a much more stable partition with periodic 
dynamics is found. This simple argument predicts that, 
under the action of noise, the system would wander be- 
tween chaotic Milnor attractors until it eventually finds a 
robust periodic attractor. If this is indeed so, chaotic at- 
tractors would always represent mere transients for suffi- 
ciently weak noises. Nonetheless, to test such hypothesis 
much longer iterations are apparently needed. 

Thus, our numerical analysis of globally coupled logis- 
tic maps has shown that the collective dynamics of this 
system in the partially ordered phase of type II is char- 
acterized by the presence of very long transients. The 
asymptotic states of the system in this parameter region 
are, however, the same as in the ordered phase and in- 
clude only a small number of synchronous clusters. This 
conclusion holds even when a small noise, eliminating 
spurious attractors, is introduced. We have also per- 
formed the analysis of dynamical transients in the in- 
termittent phase, i.e. at the interface between the or- 
dered and the turbulent phases, and have found (to be 
separately published) that in this region the coexistence 
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